1. Import data procedure
library(haven)
mining <- read_sas("http://www.principlesofeconometrics.com/sas/mining.sas7bdat", 
    NULL)
  1. Wrangle the data

library(dplyr)

mining %>%
 filter(!(YEAR %in% c("1972", "1973", "1974", "1975", "1976", "1977", "1978", 
"1979"))) %>%
 filter(QTR %in% "4") %>%
 filter(PRO >= 100.4 & PRO <= 118.6) %>%
 filter(POW >= 
    67L & POW <= 105L)
  1. ggplot2 data viz
p <- ggplot(mining) +
 aes(x = PRO, y = POW) +
 geom_point(size = 1L, colour = "#0c4c8a") +
 theme_minimal()

p

  1. Plotly

library(plotly)
package ‘plotly’ was built under R version 4.0.3Registered S3 method overwritten by 'data.table':
  method           from
  print.data.table     
Registered S3 method overwritten by 'htmlwidgets':
  method           from         
  print.htmlwidget tools:rstudio

Attaching package: ‘plotly’

The following object is masked from ‘package:ggplot2’:

    last_plot

The following object is masked from ‘package:stats’:

    filter

The following object is masked from ‘package:graphics’:

    layout
ggplotly(p)
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKMS4gSW1wb3J0IGRhdGEgcHJvY2VkdXJlCgpgYGB7cn0KbGlicmFyeShoYXZlbikKbWluaW5nIDwtIHJlYWRfc2FzKCJodHRwOi8vd3d3LnByaW5jaXBsZXNvZmVjb25vbWV0cmljcy5jb20vc2FzL21pbmluZy5zYXM3YmRhdCIsIAogICAgTlVMTCkKYGBgCgoyLiBXcmFuZ2xlIHRoZSBkYXRhCgpgYGB7cn0KCmxpYnJhcnkoZHBseXIpCgptaW5pbmcgJT4lCiBmaWx0ZXIoIShZRUFSICVpbiUgYygiMTk3MiIsICIxOTczIiwgIjE5NzQiLCAiMTk3NSIsICIxOTc2IiwgIjE5NzciLCAiMTk3OCIsIAoiMTk3OSIpKSkgJT4lCiBmaWx0ZXIoUVRSICVpbiUgIjQiKSAlPiUKIGZpbHRlcihQUk8gPj0gMTAwLjQgJiBQUk8gPD0gMTE4LjYpICU+JQogZmlsdGVyKFBPVyA+PSAKICAgIDY3TCAmIFBPVyA8PSAxMDVMKQpgYGAKCjMuIGdncGxvdDIgZGF0YSB2aXoKCmBgYHtyfQpwIDwtIGdncGxvdChtaW5pbmcpICsKIGFlcyh4ID0gUFJPLCB5ID0gUE9XKSArCiBnZW9tX3BvaW50KHNpemUgPSAxTCwgY29sb3VyID0gIiMwYzRjOGEiKSArCiB0aGVtZV9taW5pbWFsKCkKCnAKCmBgYAoKNC4gUGxvdGx5CgpgYGB7cn0KCmxpYnJhcnkocGxvdGx5KQpnZ3Bsb3RseShwKQpgYGAK